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, Abstract. We present a new modeling tool for planetary nebulae (PNe), 

C/3 ' based on 3D photoionization calculations. From models for two theoreti- 

. cal PNe, we show that the enhancement in the equatorial zone observed in 

{Sj ' several PNe is not necessarily due to a density gradient, as usually inter- 

preted. We study the morpho-kinematic properties of the PN NGC 3132 
^ ; and show that a bipolar Diabolo shape successfully reproduces the ob- 

served images, as well as the high resolution observations. 

: 

3" 1. Introduction 

ON 

I Planetary Nebulae show different morphologies, but since the works of Kwok et 

al. (1978) and Balick (1987), it is possible to reproduce their observed shape 
basically with two types of geometry: spherico-elliptical ones and bi-polar ones 
(also called butterfly). The correlation between morphological types and other 
properties linked to the central star mass set constraints on stellar evolution 
theories (e.g. Stanghellini et al. 1993, Corradi & Schwarz 1995). Here a new 
^ \ modeling tool for studying the morphology and kinematics of PNe is presented. 

This tool is first used to show that an equatorial brightness enhancement does 
not always correspond to a density enhancement. Secondly we present two 
models for NGC 3132. The first one is an ellipsoidal shell, as suggested in 
the literature, but which neither reproduces the density variation along the 
nebula, nor the velocity profiles. A Diabolo shape is successful in reproducing 
consistently all the observations. 

2. A new modeling tool 

The results of a 3D photoionization code (Gruenwald, Viegas & Broguiere 1997) 
are analyzed through an IDL (RSI) code, to generate images to be directly con- 
fronted with various kinds of observations. Rotations, projections, line imaging, 
diagnostic line ratios maps, velocity line profiles and PV-diagrams can be per- 
formed. 

3. The luminosity equatorial enhancement 

We computed two models of PNe with the same ionizing spectra and the same 
mean density, the only difference being the matter distribution: the first one 
(Model 1) have a spherical cavity surrounded by a shell with a density gradient 
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Figure 1. Emission- line imaging from Model 1 (four left panels) 
and Model 2 (four right panels). For each model, H/3, [O III]A5007, 
He IIA4686, and [N II]A6583 intensity maps are shown. 



(equatorial enhancement) while the second one (Model 2) is a constant density 
ellipsoidal shell. The global spectroscopic results are similar. The theoretical line 
images are displayed in Fig. [l], showing an equatorial luminosity enhancement. 
One way to distinguish between the two geometries is a density map determined 
e.g. by the [S II]6718/6732 ratio image. See Morisset et al. (2000) for more 
details. 

4. A 3D model for NGC 3132 

We applied 3D code to model NGC 3132, classified as an elliptical PN. We 
adopt the central star and gas properties used by Bassgen, Diesch &i Grewing 
(1990) and perform two models differing in the gas distribution. We show that 
an ellipsoidal shell sucessfully reproduces the low resolution observations (line 
images and the global spectroscopy). However, high (spatial and spectroscopic) 
resolution observations are not reproduced. Observations of [S II] 6718/6732 line 
ratio indicate a decrease of the eletronic density towards the center of the PN 
(Juguet et al. 1988). The density variation modeled with the ellipsoidal shell 
(Fig.2, dashed line) can not reproduce the observed density decrease. Choosing 
a Diabolo morphology (Fig. right side) we successfully reproduced the central 
density hole (solid line, left side). In Fig. ^ the line images of NGC 3132 from 
HST observations (top four panels) are compared to the results of our Diabolo 
model (bottom four panels). The ellipsoidal image as well as the ionization 
structure are well reproduced with a 40°rotation of the symmetry axes relative 
to the line of sight. Adopting a linear velocity law, we can model the emission 
line profiles observed through various apertures. The observed [O III] profiles 
(Sahu & Desai 1986) can be compared to the theoretical ones obtained with the 
Diabolo model using their aperture size and positions ( Fig. |j). Notice that 
the asymmetrical profiles observed at the external part of the PN (four bottom 
panels) are not obtained with the ellipsoidal shell model. See Monteiro et al. 
(2000) for more details 

5. Conclusions 

We show that the lack of high resolution observations and of a 3D modeling tool, 
can lead erroneous conclusions about the morphology of PNe. In particular, we 
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Figure 2. Left side: Variation of the electronic density along the 
PN, for the Diabolo model (solid line) and the ellipsoidal model (dashed 
line). Right side: Diabolo shape. Only high density gas (nn = 1300 
cm -3 ) is shown. The background density is 300 cm -3 . 
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Figure 4. [O III] Velocity profiles from our Diabolo model, obtained 
at the same positions as observed by Sahu & Desai (1986). 



showed that NGC 3132 was misclassified as an elliptical PN. A Diabolo shape 
is successful in reproducing: 1) the density decreasing in the central direction of 
the PN, 2) the asymmetric [O III] velocity profiles observed at the external part 
of the PN. Both observations are not reproduced by an ellipsoidal shell model. 
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